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\Q ; ABSTRACT 

£NJ . Mixed cold-hot dark matter cosmological models (CHDM) with ^Itot = lj a P _ 

^ ' proximately flat initial spectrum of adiabatic perturbations and 1, 2 or 3 types of 

on : massive neutrinos are compared and tested using recent observational data. The 

models with 2 or 3 neutrino types of equal mass permit as the best fit larger values 
^vq . of both the Hubble constant {Hq < 60 for 2 types, Hq < 65 for 3 types) and the 

total Vt y (up to 0.3 for 3 types) than the model with 1 massive type. Also, they 
. have less problems with abundances of early compact objects including Ly — a 

ON ■ clouds. 

1. Introduction 

o : 

& , If the classification of different cosmological models by a number of additional 

fundamental parameters used in them to explain all observational data EJfl is applied 
to inflationary models, then the model of the first level having only one fundamental 
parameter - an amplitude of perturbations - appears to be the CDM model with 
the approximately flat (Harrison-Zeldovich, or n ~ 1) spectrum of initial adiabatic 
perturbations. Because of theoretical considerations and observational uncertainties, 
it is better to include "weakly-tilted" models with |n — 1| < 0.1 into this class, too. 

At present, it is clear already that predictions of this model, though being not 
far from observational data (that is remarkable for a such a simple model with only 
one free parameter), still definitely do not agree with all of them. Namely, if the free 
parameter is chosen to fit the data on scales (100 — lOOO)/^ 1 Mpc, discrepancy of 
about twice in perturbation amplitude arises on scales (1 — 10)h^ 1 Mpc, and vice 
versa (h^ = Hq/50, where Ho is the Hubble constant in km/s/Mpc). Thus, models 
of the next (second) level having one more additional constant have to be considered. 
Among the best of such models is certainly the mixed cold+hot dark matter model 
(CHDM) i, for recent analysis see our previous papers (hereafter PS1 § and PS2 I), 
as well as tffl and references therein. In this model, the hot component is assumed to 
be the most massive of 3 neutrino species (presumably, r-neutrino) with the standard 
concentration following from the textbook Big Bang theory, while masses of the other 
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two types of neutrinos are supposed to be much less and, therefore, unimportant for 
cosmology. Then the only new fundamental parameter is the neutrino mass m v . 

Still the CHDM model with n ~ 1 is not without difficulties. The main of them 
is connected with later galaxy and quasar formation in this model as compared to 
the SCDM model. As a result, only a small region in the Hp — Q u plane remains 
permitted (PS1 a, some authors have a more pessimistic view 0). The analysis of the 
possibility of non-flat (but power law) initial spectra (the article Q for n < 1, PS2 □ 
for detailed analysis of n > 1 case) leads to the conclusion that the CHDM works the 
best with nearly flat n « (0.95 — 0.97) initial spectra, which follow from the simplest 
inflationary models. Therefore, allowing for this degree of freedom does not improve 
the fit to observational data. 

In this paper we consider how the cosmological predictions of the CHDM change if 
not only one but two or even all three types of neutrino are massive and contribute to 
the present density of hot matter. The total energy density (in terms of critical one) 
of hot component in this general case is then ^m Vj = 23AVL v h\ eV for T 7 = 2.73K. 

This version of CHDM, of course, also corresponds to the next level of complexity 
in our classification until the mass ratio of different neutrino will be either confirmed 
in laboratory experiments, or theoretically derived from some underlying theory. We 
shall not cover all possible combination of the neutrino masses, noticing that the case 
of three equal masses m Vl = m V2 = m„ 3 is the one mostly different from a standard 
model with one massive neutrino. The predictions for all other sets of m Vi (with the 
same total Vt v ) lie "in between" these two models. There is no experimental evidence 
for all three types of neutrino to have comparable mass. Thus the model with m ui = 
m U2 = m„ 3 is mostly interesting as a limiting case. More realistic is the assumption of 
two type of neutrino to have mass in cosmologically interesting range while the rest 
one (electron neutrino) has a much smaller mass. Therefore we restrict ourselves with 
two models m Vl = 0, m V2 = m U3 and m Vl = m V2 = m V3 and compare them with the 
original CHDM with one massive neutrino. We denote by N v the number of (equally) 
massive types of neutrino which mass is now m v = 23AVt v / N u h\ Q eV. 

2. CHDM Linear Perturbation Spectrum 

The transfer function C(k) for the CHDM with one massive neutrino (N u = 1) was 
determined in our previous paper □ numerically by solving a system of the Einstein- 
Vlasov equations for the evolution of adiabatic perturbations with a neutrino com- 
ponent treated kinetically and a CDM component as dust. The correspondent fitting 
formula is given in □. Results of similar computations for the case N v = 2,3 are given 
in Fig.l in comparison with N v = 1 curve for a specific choice fl u = 0.2. 

The main difference between models comes from the fact that the same total 
mass density fraction Q u is achieved with different masses of neutrino. During the 
matter dominated stage the presence of hot component distributed at the begin- 
ning of this stage uniformly due to free-streaming lead to the retarded growth of 
perturbation in total density on modes kR nr < 1. Here R nr is the comoving size 



Figure 1: Upper-left panel: C(k) transfer function for Q u = 0.2 and N v = 1,2,3 massive 
neutrino types. Other panels: Restrictions in the Hq — VL V parameter plane for N u = 1 (upper 
right), N v = 2 (lower right) and N v = 3 (lower left panel), following from: a) fit to Stromlo-APM 
counts-in-cells. Solid lines correspond to % 2 = 2, 7 contours; b) the <Jg condition. The values 
cr 8 < 1, 0.67 for Q rm s-ps = 14.3 fiK are achieved left to the dashed (correspondingly rightmost 
and leftmost) lines. The middle dashed line corresponds to <Jg = 1 if Qrms-PS = 19 /iK; c) the 
combination of the <Tg < 0.67 condition with quasar and galaxy formation conditions (dotted 
lines). The allowed region lies below the upper dotted line if the limitation set for objects of the 
mass M = 10 11 Mm is used and below the lower one if M = 10 12 M ra . 



of horizon at the moment when neutrino become nonrelativistic. For the mode k 
this period continues until k becomes smaller that the effective neutrino Jeans scale 
fc 2 = 4:7rGp tot /(dp u /dp u ) oc m 2 t~ 2 / 3 when perturbation in neutrino density reach the 
magnitude of inhomogeneities in cold component (see discussion in PS1 a). This pe- 
riod is larger for smaller mass, and as the result, the amplitude of perturbation at 
the present moment is decreased even further from SCDM value. This explains why 
the transfer function Cchdm is ~ 20% lower for larger N u at the intermediate scales 
R~l < k < fcj(io) ~ H^i/^io Mpc -1 . For wavelength shorter that the present Jeans 
scale k > kj(t ) the perturbation in neutrino density did not start to grow up to the 
present moment and no additional dependence on m u appear. In this limit, as Fig.l 
demonstrates, the resulting amplitude is approximately the same as far as the total 
amount of neutrino Q u is the same. 



3. Confrontation with Observational Tests 



In Fig.l we present the restrictions in the Vt v — H parameter plane which follow 
from several observational tests shown to be the most illustrative in our previous 
papers. These tests ) value of the total rms mass fluctuation ag at R = 

16/i5 1 Mpc (er 8 < 0.67 based on cluster abundance data @) as follows from COBE 
measurement of DMR anisotropy (AT/T) which we adopt as Q r ms-ps — (17.4 ±3.1) 
/iK El; b) The Stromlo-APM counts in cells £-1 which limits the slope of spectrum 
over the range (20 - 150) /igo 1 Mpc; c) density of quasars at high redshifts z ~ 4 
Necessity to produce sufficient number of quasars in the model puts a lower limit on 
mass fluctuations at oc IMpc scale. We follow @ where the quasar density criterion at 
z = 4 was formulated as a restriction on a mass fraction in bound objects with mass 
larger 10 11 M Q : /(> 10 11 M Q ) > 10~ 4 . Similar estimate for a fraction of mass in 
large galaxies gives /(> 10 12 M Q ) > 10 -5 at the same redshift z = 4. We put greater 
weight on this second, more restrictive limit, because, remarkably, it gives similar 
restrictions on the CHDM model as the analysis of the abundance of damped Ly — a 
absorption systems in high-redshift quasar spectra 00. 

The upper-right panel of Fig.l presents the fit, — Hq parameter plane with obser- 
vational restrictions for standard CHDM with n = 0.95. This plot was discussed in 
detail both in PS1 a and PS2 H, showing that CHDM model parameters are restricted 
to the narrow range of a low Hubble constant H < 55 km/s/Mpc and the neutrino 



fraction Vt u = 0.17 — 0.28 for Hq = 50. The first of these limits reflects a problem 
with unavoidable high mass fluctuations at the 16/i^q 1 Mpc scale (note that to set this 
limit we use a lower bound on (AT/T), Q rms -ps — 14.3 fiK although recent analysis 
suggest higher value Q r ms-ps ~ 20 /iK), as well as a wrong shape of the perturbation 
spectrum over the I = (20 - 150) Mpc interval if H is high. The upper bound on 
Q u comes from the combined quasar and as conditions that implies that the slope of 
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the spectrum in the scale range (0.7h 5Q — 16h§ ) Mpc cannot be too steep. 

The bottom panels of Fig.l correspond to two- and three- massive neutrino models 
N„ = 2, 3. The additional decrease of the transfer function in the range O.OlMpc -1 < 
k < IMpc -1 while keeping both larger and smaller scales unchanged allow us to 
noticeably relax the restrictions on the CHDM models. The major change is the 
possibility to accommodate the higher values of Hubble constant. Now both the 
shape of spectrum better fits count-in-cells and the as limit, understandably, becomes 
less severe. For an extreme example N v = 3 we may have Hq as high as 70 — 75 
km/s/Mpc. For N u = 2 the tests can be satisfied if Hq < 65 km/s/Mpc, although 
more comfortable fit to data is achieved for H < 60 (N v = 2) and Hq < 65 (N u = 3). 
Let us stress, however, that CHDM by itself does not address the problem of the age of 
the Universe, which remains the major theoretical objection to Hq > 65 km/s/Mpc. 

Additionally a somewhat higher values of Q u are now allowed. Although for 
Hq — 50 the boundaries on total neutrino fraction remain unchanged Vt v = 0.15—0.25, 
for H = 60 and N v = 2, 3 we have from Fig.l 0.2 < Q u < 0.3. This change is rather 
significant as a limit on the sum of neutrino masses: 3.5eV < m Vi < 6eV if we take 
Hq = 50 and 6.7eV < ^m u . < lOeV for Hq = 60. On the other hand, the mass of a 
single neutrino type is approximately the same m Ui ~ 2 — 5 eV. 

We can conclude that having several types of neutrino with cosmologically signif- 
icant mass of few electron-volts expands the boundary of parameter space for CHDM 
model. The allowed region of Q u — Hq plane now include not only the small area 
of Hq ps 50, VL U ~ 0.2 but extends to H ps 60, Q u ~ 0.25 for the most physically 
interesting version with two massive neutrino, and as far as Hq k, 70, Vt v ps 0.3 for 
the extreme case of three equally massive types. Similar results for the N v = 2 case 
were recently presented in EJ. 
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